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ABSTRACT 

NGC 6240 is a rare object in the local universe: an active merger remnant viewed at the point of 
merging where two active galactic nuclei are visible. We present IRAC data of this object, providing 
high sensitivity maps of the stellar and PAH distribution in this complicated system. We use pho- 
tometry to analyze the variation in these distributions with radius and provide an SED in the four 
IRAC bands: 3.6 /im, 4.5 /im, 5.8 /im and 8 /im. We fit the radial profiles of the 3.6 /im band to r * and 
exponential profiles to evaluate the structure of the remnant. Finally, we compare the IRAC images 
with multi- wavelength data and examine how outflows in the X-ray, Ha and CO correlate with 8 fim 
emission. The results support the general picture of NGC 6240 as a system experiencing a major 
merger and transitioning from a disk galaxy to a spheroid. The sensitivity of IRAC to low surface 
brightness mid-infrared features provides detailed information on the extended distributions of stars 
and dust in this rare system. 

Subject headings: galaxies: photometry, galaxies: active, galaxies: irregular, infrared: galaxies, galax- 
ies: interactions 



1. INTRODUCTION 

Due to its proxim ity, NGC 6240 (z = .0243 
iSolomon et all |l99l . L> = 98 Mpc for Ho = 75 
km s~^ Mpc~^) is one of the brightest objects in the 
IRAS All-Sky Survey and was one of the early iden- 
tified LIRGs with a total i nfrared luminosity of 7.1 x 
10"io (Wright et al.llT98l . Optical studies of NGC 
6240 reveal a morphologically irregular galaxy with 
clear tidal extensions, two nuclei and lar ge dust lanes. 



stron g ly suggesting a mer g ing systeni (IZwickv et al 
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119831: iKe el" 1990; Gcrssc n et al. 
20041) . Near-IR data (iThronson et al. 1990: Dovon et al 
1991 iScoville et all I2OOOI: iMax et all I2005D an d ra 



dio continuum (ICarral et al.l 1990l: Bales et al.l 119901 : 
IColbert et al.l 119941 : iBeswick et al.ll200lD images confirm 
the double nuclei and show that they have consider- 
able substructure. The observed distance between the 
the nuclei depends on the wavelength at which they 
are observed, indicating large amou nts of dust extinc- 
tion in the center of the galaxy (jSchulz et al.l 119931 : 
iGerssen et al.l [2004) . Ha data shows stro ng evidence 
for n uclear winds, outflows and bubbles (Gerssen et al.l 
|2004[ ). X-ray studies show a very active galaxy with 
two nuclei dis playing AGN characteristics , indicating two 
black holes (llwasawa fc Comastril 119981 : IVignati et al.l 
Il999t iKomossa et al.ll2003D . Sub- millimeter array data 
reveal a molecular gas disk of mass « 10^ Mq in be- 
tween the two nuclei with the possibility of outflow s 
(|Wang et al.lll991l: iTacconi et allll999l : llono et al.|[2007l ) . 
Merger remnants with double nuclei are rare in the 
nearby universe, suggesting that we are catching this ob- 
ject in a very short lived phase of its evolution. In addi- 
tion, its infrared luminosity places it on border between 
the LIRG and a ULIRG classes. Though classified by 
their luminosity, the LIRG and ULIRG classes show dif- 
fering characteristics besides their luminosity, including 
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a higher occurrence of AGN and dist urbed mo rphologies 
in ULIRGs (Veilleux et al. 2002; Sand ers fc Ish ida 2004). 
Given that NGC 6240 already shows many of these fea- 
tures, we may be observing a transition object between 
the LIRG and ULIRG classes. 

Objects such as NGC 6240 are believed to be under- 
going a merger between two comparable mass galaxies 
which is driving the remnant through evolution from 
a Seyfert type active galaxy to an AGN and finally 
a quasar (iNorman fc Scovilld[l988l : ISanders et al.lll98a 
iHopkins et al.l l2006[ ). Models demonstrate a cycle in 
which a merger drives gas inflows that create starbursts 
and feed central super-massive black holes, triggering an 
active stage for the black hole. As the AGN "turns on" 
and contributes significant radiative feedback, it clears 
out the gas around the nuclei, and the AGN becomes vis- 
ible as a quasar. However, this also depletes its fuel, leav- 
ing the AGN to die back to a quiescent phase. The cycle 
may start ove r again when the gala xy goes through an- 
other merger. iHopkins et al.l (|2008al la) show this cycle is 
consistent with many observed properties of quasars and 
elliptical galaxies. This evolutionary sequence is short 
{r-^ few Gyr) compared to the lifetime of a spiral or ellip- 
tical galaxy, so there are relatively few of them to study, 
and even fewer at the particular stage of merging where 
two AGN are observable. This makes them particularly 
interesting objects to stu dy in detail. 

We add Spitzer/IRAC (|Fazio et al.ll200l data of NGC 
6240 to the existing observations of NGC 6240. Most in- 
frared studies of NGC 6240 have focused on the nuclei 
(iKlaas et all 120011 : iLutz et al.ll2003l : lArmus et"al] 120061 : 
Egami et alJI . 20061 1. We draw on IRAC's high sensitiv- 
ity and multiple mid-infrared bands to investigate how 
stars and dust are distributed throughout the remnant. 
The IRAC band-passes single out emission from stars 
(3.6 /im) and polycyclic aromatic hydrocarbon emission 
(PAH, a component of dust, 8/im), allowing us to dissect 
the structure of NGC 6240 by examining the distribution 
of stars and ISM separately (Pahrc et al. 2004). All four 
IRAC bands are imaged on the same scale so that these 



components can be easily compared. We derive light pro- 
files in each of the IRAC bands and calculate the SED. 
We also compare IRAC data to archive optical, X-ray 
and radio data. In fj2] we describe our observations and 
in Sj3]we describe the IRAC images and how they com- 
pare to data from other wavelengths. In f|4] we present 
our photometry of the IRAC data and in [JS] we make 
some concluding remarks. 

2. OBSERVATIONS 



Ob servations were carried out with IRAC (^ Fazio et al.l 
|2004| ) as part of the Spitzer GTO program on 13 August 
2004 in the 3.6, 4.5, 5.8 and 8.0 /im bands. Because of 
NGC 6240's high infrared luminosity, the data was taken 
in high dynamic mode. Sixty dithered images with an ex- 
posure time of 0.6 s were taken in each filter in addition 
to 120 dithered images with an exposure time of 10.4 s. 
The short exposure images allow checks against satura- 
tion and bright source effects. The short and long ex- 
posures were mosaic-ed separately giving total exposure 
times of 10.8 mins and 36 s in the long and short expo- 
sure mosaics, respectively. The IRAC post pipeline Ba- 
sic Calibrated Data (BCDs, pipeline version 14.0.0) were 
corrected for artifacts with an artifact mitigation code 
(ICarevll2007[) a r id wer e mosaiced with IRACproc (vl4.1, 
ISchuster et all (|2006[ )). The pixel size of the BCDs is 
1.22". The pixel size of the mosaiced images is 0.86", 
this is the minimum re-sampling needed to ensure that 
the IRAC PRF is not under-sampled even at 3.6 and 4.5 
/xms. The FWHM of the IRAC point response functions 
are all ~ 2.0", which corresponds to 950 pc at the dis- 
tance of NGC 624 0. The double AGN of NGC 6240 are 
~ 1 kpc apart (e.g. lGerssen et alll2004f) . so we cannot re- 
solve them and will concentrate on the extended emission 
of the remnant. 

We compare the IRAC data to archival data from many 
sources including Chandra Soft X-ray, Hubble Space 
Telescope Advanced Camera for Surveys (ACS), Very 
Large Array (VLA) radio continuum and sub-millimeter 
array (SMA) CO 3-2. The Chandra observations of 
NGC 6240 were carried out over 37 ks on 2001 July 
29 (OBSID 1590) using the ACIS-S instrument. The 
pixel size of the ACIS detector is 0.49". We analyzed 
the data using CIAO^ data reduction and analysis rou- 
tines. To detect the low contrast diffuse emission sur- 
rounding the bright nuclear region we ap plied the CIAO 
adaptive smoothing tool csmooth ( Ebcli ng et al.l I2006D . 
Csmooth preserves the mult i-scale spatial signatures ari d 
the associated counts fe.g. iKarovska et 311120021 . l2007f ). 
The smoothing is achieved by convolution with a Gaus- 
sian kernel which increases from a small initial smoothing 
scale until the total number of counts under the kernel 
exceeds a value that is determined from a preset signifi- 
cance and the expected number of background counts in 
the kernel area. The maximum smoothing scale in the 
image is 5" . 

Hubble Space Telescope ACS data in the B broadband 
filter and Wide Fiel d Planetary Camera (WFPC2) Ha 
(|Gerssen et al.ll2003 ) data were obtained from the NASA 
Multi-mission Archive at STScI (MAST)3. VLA radio 

^ CIAO is the Chandra In teractive Analysis of Observation's 
data analysis system package l |http://cxc. harvard. edu/ciaO| ) 
^ http://archive.stsci.edu 



continuum data centered at 21-cm was obtained from 
the N ASA extragalactic database (NED'*: ICondon et al.l 
|1996[ ) and sub-millimeter array data of CO 3 -2 was ob- 
taine d with the Submillimeter Array (SMA, llono et aLJ 
l200l . 

3. RESULTS 
3.1. IRAC Imaging 

In Figured) we show the IRAC images in each bandpass 
along with the ACS B-band, ACS I-band and Chandra 
soft X-ray images. Contours are overlaid on the IRAC 
images to highlight structure. The lowest contour is set 
at 4cr above the background, which is .16, .33, 5.06, 9.68 
MJy/str in the 3.6, 4.5, 5.8 and 8.0 /im bands respec- 
tively. Four features extending off the remnant are la- 
beled in the 3.6 /xm panel of Figure [Tj They point ap- 
proximately south, north, southeast (then curves around 
to the south) and northeast (smaller) and are labeled 
as such. The 8.0 /im image is dominated by "spikes" of 
emission to the east and west of the galaxy. These are due 
to the bright central source illuminating the IRAC 8.0 ^m 
point response function and are not true features. They 
are clearly seen in the PRF of the 8.0 /im band shown in 
Figure[2l panel (b). There is also some contribution from 
the PRF to the 8 ^m image along the north-south axis of 
the remnant. Three low intensity features that are not 
due to the PRF are indicated with arrows on the 8/im 
panel in Figure [T] a dust compliment to the southeast- 
ern feature, a bright source on the southwest side of the 
remnant and the dust lane to the north. Finally, note 
that the bright source to the northeast of the remnant is 
a foreground star. 

In an attempt to remove the PRF artifacts and clar- 
ify the dust features in the 8 /im image, we subtract off 
the PRF of the central emission, treating it as a point 
source (Figure Hj) . The magnitude of the subtracted ob- 
ject is ~ 6 (Vega). A few artifacts are left in the center 
after the subtraction, and these pixels are set to an av- 
erage value of the pixels around them. In addition, the 
reflection artifacts shown on cither side of the remnant 
are cleaned by averaging over columns of bright and dark 
pixels. This means that for these six pixels on either side 
of the remnant, features are smoothed in the NE-SW 
direction. The final image is shown in Figure [2]i. Fea- 
tures far from the central source of the remnant should 
be largely unaffected by remaining artifacts, but any re- 
maining axisymmetric features should not be trusted. 
This image is never used for photometry, but it allows 
us to exploit IRAC's strengths to examine the morphol- 
ogy of the extended emission in NGC 6240. Figure [2] 
shows the original 8/im image, the 8/im PRF and the 
final image. Dust features along the northern, southern 
and southeastern extensions are now clearly visible and 
show noticeable substructure. The northern extension in 
particular appears very clumpy. The southeastern exten- 
sion is pronounced in 8 /im and seems to be a long, thin 
tidal tail. The southern extension is quite faint compared 
to the other features in the 8 /xm and curves such that the 
tip is pointing to the southeast. Along the southeastern 
side of the remnant, between the southern and southeast- 
ern extensions, dust filaments are present, suggestive of 
entrained dust being blown out of the remnant. 

^ http://ncdwww.ipac.caltech.edu/ 




Fig. 1. — Multi-wavelength data of NGC 6240. On top are the four IRAC broadband images with contours. North is up and east is to 
the left. In all images the contours are logarithmically spaced with the lowest contour at Acr above the background. The minimum contours 
are .16, .33, 5.06, 9.68 MJy/str and maximum contours are 68.79, 88.59, 249.2, 400.1 MJy/str for the four bands respectively. In the 5.8 /^m 
images additional contours are added between the first two contours at 5.5, 6, and 6.5 MJy/str. In the 8.0 /im image additional contours 
are added between the first two contours at 10, 11 and 12 MJy/str. A 15"scale is shown in the image. At the distance of NGC 6240, this 
corresponds to 7.1 kpc. Note that the 5.8 and 8.0 /^m bands show prominent PRF effects, see the PRE in Figure [2] On the bottom are 
ACS B, I and Chandra soft X-ray images of NGC 6240 at the same scale. Logarithmic contours are shown on the X-ray image. 

bands, therefore these bands trace emission from com- 
plex dust molecules. The standard IRAC coloring is 
used in the three-color image (3.6 /im - blue, 4.5 /im - 
green and 8.0 /im - red) so that stellar populations ap- 
pear blue, PAH emission appears red and combinations 
of the two appear white. The most striking aspects of the 
IRAC 3-color image are the differences across the north- 
ern and southeastern features - highlighting segregation 
of the stellar and dust distribution. The northern exten- 
sion shows clear separation of the stars and dust - stars 
dominating on the eastern side and the dust dominating 
on the western side. The southeastern extension present 
in the optical is clear in the 3.6 /im, and the 8 /im reveals 
a dust counterpart lining the eastern side. The long, thin, 
arched appearance of this extension is characteristic of a 



^#44 



Fig. 2. — This figure illustrates how we processed the 8 fim image 
to remove the central PRF. (a) The original 8 fiva image, (b) The 
point response function of the 8.0 /im band showing the PRF spikes 
that afflict the 8 fim image. This PRF is saturated to show the 
outer PRF features, so these effects will be significantly fainter 
in the image of NGC 6240. (c) The 8 fim image with the PRF 
subtracted out. Because this is not a point source, the subtraction 
cannot be done perfectly, and some artifacts are left at the center 
of the image. The lines at the side are a well known IRAC artifact, 
(d) Finally, these artifacts are cleaned to make the image more 
visually pleasing. No photometry is done on this image, but it is 
used for creating all the following images. 

To examine how components of the remnant are dis- 
tributed, the IRAC 3.6, 4.5 and 8/im data are combined 
to create a three-color image of NGC 6240 in Figure [3] 
(left). The 3.6 and 4.5 /im bands, with less than one 
tenth the dust obscuration of the optical, trace old stars, 
despite significant dust obscuration. The 7.7 /im and 
8.6 /im PAH features are located in the 5.8 and 8.0 /im 



tidal tails. A gas/star offsets in ti dal tails are commonly 
observed in merger r emnants (e.g. iHibbard fc Yuiilll999l : 
iHibbard et al.|[2005f) and can be explained by the coUi- 
sional nature of the gas (jMihosI 12001). The right panel 
of Figure [3] is a three-color combination of the ACS B 
(blue), I (green) and PRF subtracted 8/xm (red) images. 
Note that the resolution of the ACS data is an order of 
magnitude better than the IRAC data. The yellow colors 
in the right panel correspond to areas dominated by I and 
8 /im, while the red colors indicate only 8 /im emission is 
present. Young, unobscured stellar populations appear 
blue. Again, areas with both dust and stellar emission 
appear white. Since the B and I band are subject to 



more dust obscuration than the 3.6 and 4.5 /im bands, 
only reasonably unobscured stellar populations appear 
white in the right image, while stellar populations that 
suffer more dust obscuration still appear white in the left 
image. 

The 8 /xm emission and optical dust lane correlate 
amazingly well along the northern extension, down to 
individual clumps in the 8 fj,m emission. However, the 
8 fj,m also reveals dust we would not have seen in the 
optical, such as the dust/gas companion to the SE tail. 

More subtle features are revealed in the three-color im- 
ages (Figure [3]) as well. A faint southern dust exten- 
sion, revealed by the 8/im, is spatially coincident with 
the southern stellar extension in the optical and 3.6 ^m, 
rather than being separated as in the northern and south- 
eastern extensions. This extension is very irregular, fil- 
amentary and truncated. Filamentary 8/zni emission is 
also present along the southeast side of the remnant, be- 
tween the southeastern and southern extensions. This is 
most clear in Figure [21 panel d. Nothing like this emis- 
sion is seen in the shorter IRAC wavelengths, indicating 
this filamentary dust is unacc ompanied by stars. Sim i- 
lar filaments are seen in M 82 (jEngelbracht et al] 120061 ). 
The northeastern extension also seems to be dominated 
by dust, especially along its southern edge. 

Overall, the stellar distribution is much smoother than 
the dust distribution. This may be a combination of two 
effects: differing distributions of gas and stars in the pre- 
merger galaxies and the fact that physical processes in 
mergers affect stars and gas differently. Notice, in Fig- 
ure [U that the 3.6 /im contours are evenly spaced and 
straight on the east and west side of the body of the 
remnant, while the dust obscuration of optical light cre- 
ates the irregular "bow-tie" appearance of this remnant. 
The 8 /im emission shows particular elongation along the 
southwestern part of the bow-tie, rather than following 
the smooth distribution of the stars. Each of the ex- 
tensions is more clumpy in dust emission than in stellar 
emission. Along the northern extension and the body 
of the remnant, where the morphology of the remnant 
looks like an edge-on disk, the separation of stars and 
dust could just be a geometric effect (edge-on disks show 
a prominent dust lane running along their major axis due 
to dust obscuration which appears off center if the disk 
is tilted towards or away from the observer). We use 
photometry to test this in Section |4l However, in many 
parts of the remnant, such as the filamentary emission in 
the southeastern part of the remnant, the different dis- 
tributions are likely to reflect the fundamental fact that 
mergers affect gas differently than stars, due to its coUi- 
sional nature. The clumpy and uneven dust distribution 
is a clear prediction of starburst and AGN feedback (e.g 
iSpringel et an[200l . 

3.2. Comparison to Multi- Wavelength Data 

We compile optical, soft X-ray, Ha, CO and radio data 
of NGC 6240 for comparison with the IRAC data. Fig- 
ure [4] shows the contours from the smoothed soft X-ray 
image (0.5-1.5 keV) overlaid on the IRAC data where 
the extended diffuse emission was the most prominent 
(~ 1000 total counts). Soft X-rays trace hot gas (> 10^ 
K) . The X-ray data shows an X-like structure in the dif- 
fuse emission extending to ~1 arcminute from the cen- 
tral region with an axis of symmetry along the north- 



south features of NGC 6240. The northeast-southwest 
line of the X is especially prominent. This pattern is 
probably associated with stellar outflows generated by 
the starburst in NGC 6240; soft X-ray emission from 
similar stellar outflows in s tarbursts is w ell documented 
in cases like M 82 (Grifliths et al.l 120001 ). Some of the 
emission may b e scattered X-ra.y ligh t, as in the case 
of Centaurus-A (jKarovska et al.l [20021 ). In the inner re- 
gions the X-ray shows a "cloverleaf" pattern with three 
clear bubbles, most likely associa ted with outflows gen- 
erated by the AGN ([Komossa et al. 20031 . The Ha data 
follows the inner X-ray contou rs closely and appears fil- 
amentary (jGerssen et al.|[2C)04D . another good indication 
of outflows in this remnant. 

CO (tracing molecular gas) and radio continuum emis- 
sion (roughly tracing synchrotron emission) contours are 
shown in Figure [5| The 3a contour for the Ha emission 
is shown here in green and CO data is shown in cyan . 
We only show CO emission peaks that llono et all (|2007f ) 
consider significant. Radio emission is concentrated in 
the center, but also shows a peak arcing to the west. 
All three forms of gas emission show a peak over the 
central source, but the CO and radio data also show a 
peak inside a bubble outlined by the west two filaments 
of Ha emission. This bubble is also coincident with the 
west "leaf" of the "cloverleaf" of X-ray emission shown 
in Figure [H ([Komossa et al.l [20031 ) . This has been inter- 
preted as a large outflow (hereafter referred to as the 
western outflow), wi th evidence of the colder molecular 
gas being entrained ([Komossa et al.[[2003[ : [Gerssen et al.l 
2004; lono et al. 20(3^^ 

While the stellar component traced by IRAC is col- 
lisonless, and will not be affected by these outflows, the 
dust may be affected. The 8 /im emission clearly peaks 
over the central source, but does not seem to show a 
corresponding peak over the western outflow seen in the 
Ha, CO and radio. PAH emission in the western out- 
flow region seems smooth, with no evidence for a peak 
or blowout. However, the PRF of the central source and 
effects from the PRF subtraction could affect our ability 
to see a peak or minimum in this region. The filamen- 
tary emission along the southeastern side and the patchy, 
almost filamentary emission coincident with stars in the 
southern extension are tantalizing indications of dust en- 
trainment in outflows. It is possible that the tempera- 
ture and pressures in the western outflow are so high, as 
shown by the X-ray emission there, that complex PAH 
molecules have been destroyed, while they survived in 
other, less hot outflows in NGC 6240. However, the CO 
peak in the western outflow shows some cold gas sur- 
vives here, which we would expect to contain complex 
molecules such as PAHs. Other bubbles in the X-ray and 
Ha seem to correspond well with the PAH emission. The 
filamentary emission in the southeastern part of the rem- 
nant seems to follow well from the southern Ha fllaments 
and southern "leaf" of the X-ray "cloverleaf" . The strong 
PAH emission in the northeastern extension is coincident 
with the northeast "leaf" of the X-ray "cloverleaf" and 
the northeast filaments of the Ha emission. 

4. ANALYSIS AND DISCUSSION 

4.1. Elliptical Photometry 

To quantify the distribution of stars and dust in NGC 
6240, IRAC colors were determined by fitting elliptical 




Fig. 3. — Optical and IRAC three-color images. Left: IRAC 3.6 /im in blue, IRAC 4.5 /im in green and PRF subtracted 8 fim in red. 
Right: ACS B-band in blue, ACS I-Band in green and PRF subtracted 8/^m in red. Notice how 8/im emission follows the optical dust 
obscuration very closely. 




Fig. 4. — X-ray contours (yellow) overlaid on IRAC three-color 
image. The first contour is at the 3a level, and they continue 
logarithmically. Lines are drawn to draw attention to the "X" 
pattern in the X-ray emission. 



isophotes to the remnant and examining the enclosed 
light to investigate changes in color with radius. Using 
the task "ellipse" in IRAF, elliptical isophotes were fit to 
the 3.6 fiia data and used as apertures for analysis of the 
4.5, 5.8 and 8/im bands. Removal of the 8/im PRF ar- 




FlG. 5. — The inner regions of NGC 6240 with contours showing 
different elements of the interstellar medium overlaid on the optical 
three-color image from Figure O The 3cr Ho contour is shown in 
green, while CO contours are shown in cyan and radio continuum 
contours are shown in yellow. We only plot CO peaks considered 
significant bv llono et al.l II2007I ). 



tifacts quantitatively is difficult due to the fact that the 
bright central source of NGC 6240 is not a point source, 
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so photometry was performed on the original image. To 
compare at identical angular resolu tion, each of the im- 
ages was convolved with a kernel ([Gordon et al.l l2007f ) 
which matched their PRF to the 8.0 /zm PRF before the 
photometry was applied. This corrects the image for the 
dominant central effects (radius less than « 4 kpc). The 
outer PRF "spikes" contribute an error of ^2% to the 
enclosed flux in the outer aperatures. 

The semi major axis of the innermost ellipse was 0.94 
kpc {^ 2") and it was increased by 0.94 kpc for each 
larger ellipse, which corresponds to the FWHM of the 
point response function of the 8 /xm band. For each el- 
lipse, the ellipticity, position angle and central coordi- 
nates of the ellipses were allowed to vary at the prede- 
termined semi-major axis. The ellipticity of the ellipses 
increases slowly from an initial value of .25 to .5 at 8 
kpc, then remains approximately constant. The posi- 
tion angle of the ellipses varies between 80 and 90 de- 
grees, so the semi-major axis is approximately along the 
remnant's north-south extensions. The maximum radius 
was the largest radius ellipse fit before the neighboring 
foreground star was enclosed by the ellipse: 20.75 kpc 
(44"). Despite NGC 6240 's irregular morphology in the 
optical, fairly concentric ellipses were fit, reflecting the 
smoother emission in the 3.6 /im band than in the optical. 
These apertures were then used to evaluate the bright- 
ness as a function of radius in all four bands. A constant 
background was subtracted from all channels. The back- 
ground was calculated by taking a 115 kpc by 60 kpc 
box around the remnant and running 3cr clipping itera- 
tions until the mean background value converged. The 
dominant error is in the ellipticity of the fitted ellipses. 
Uncertainties in the fitted ellipses were determined by 
the "ellipse" task from the uncertainties in the harmonic 
fits and final error bars were determined by changing the 
ellipticity of the ellipse by this error and measuring the 
change in the intensity enclosed by the ellipse. Statisti- 
cal errors due to pixel noise are negligible. Errors due to 
the calibration of the original BCDs are discussed in the 
IRAC data handbook. As a guide, IRAC photometry is 
only accurate to a few percent due to calibration errors, 
though this will only affect absolute fluxes. Finally, to 
check for saturation in the 8 /zm band, we performed the 
same analysis on the short exposure mosaic and found 
virtually identical results in the inner regions. 

Using the largest elliptical aperture described we mea- 
sured the total fiux of NGC 6240 to be .10, .09, .21 and 
.49 Jy in the 3.6, 4.5, 5.8 and 8/im bands respectively. 
This SED is shown in Figure [6] with the Spitzer/Infrared 
Spectrograph (IRS) data for NGC 6240 obtained by 
lArmus et al.l |2006). Though this IRS data only covers 
wavelengths longer than 5.0 /im, it agrees with the IRAC 
SED to within .1 Jy at these wavelengths. The IRAC 
SED also ag rees well w i th th e integrated ISO spectra 
presented bv lLutz et al.l (|2003f ). Note that bands 3 and 
4 show strong PAH features, indicating the overall SED 
is weighted strongly towards the central dust-dominated 
regions of the remnant. 

The flux as a function of radius for each band is shown 
in Figure [T] The total flux in each band slowly levels 
off at large radii, indicating that the largest aperture en- 
closes most of the flux from the galaxy, and each curve 
is normalized to the total flux in that band. The longer 
wavelength bands gain flux much more rapidly in the cen- 
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Fig. 6. — The integrated spectral energy distribution of NGC 
6240 in the four IRAC bands. For comparison, IRS data is plotted 
for wavelengths longer than 5.0 ^m (Armus ct al. 2006). The IRS 
data is not scaled. The IRAC bandpasscs are also plotted as the 
colored lines (red - 3.6, green - 4.5, yellow - 5.8 red - 8.0 ^m) 
where the scale is correct but the units are number of electrons per 
incoming photon. 
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Fig. 7. — Enclosed flux at a given radius derived from fitted 
isophotal ellipses for NGC 6240: blue - 3.6 fim, green - 4.5 ^m, 
yellow - 5.8 fim and red - 8.0 fj,m. Curves are normalized to the total 
flux of the remnant in each band. Note that at longer wavelengths, 
the emission is more concentrated. 

tral few kpc than the shorter wavelength bands, which 
indicates that the flux from dust emission in the galaxy is 
more centrally concentrated than that from stellar emis- 
sion. A good way of examining the changing colors in the 
remnant is to calculate the ratios of the 5.8 to 3.6 /im col- 
ors and the 8.0 to 4.5 /im colors, effectively calculating 
at how red or blue the remnant is in the mid-IR colors 
as a function of aperture radius. The colors calculated 
from the innermost aperture of our photometry, enclosing 
both nuclei of NGC 6240, are the reddest in 5.8/3.6 /im of 
any of the apertures and the 5.8/3.6 /im emission always 
decreases with increasing aperture radius. In 8.0/5.8 /^m 
emission the colors increase to a "knee" at 2.8 kpc, where 
the aperture encloses the majority of the 8 fim emission. 



indicating the large concentration of dust emission in 
the center of the remnant. Then the colors decrease in 
8. 0/4.5 /xm wit h large r aperture. 

ISaiina et al.l ()2005f ) use SED modeling to determine 
the color ratios of galaxies dominated by continuum, 
PAH and stellar sources. They find that objects domi- 
nated by different sources of emission lie in well differen- 
tiated parts of the plot, supporting observatio nal results 
(|Lacv et al.l 12004 iHatziminaoglou et al.ll2005t ), but that 
the populations overlap at their edges. They caution that 
determining the characteristics of a galaxy using this di- 
agram is not necessarily appropriate for samples smaller 
than 10, due to the risk of the galaxy being an outlier, so 
drawing strong conclusions from NGC 6240's location in 
this plot would be unwise. However, in their modeling, 
our IRAC colors for NGC 6240 place it on the border 
between the PAH dominated and continuum dominated 
regimes, leaving the possibility that continuum emission 
from the AGN in NGC 6240 is seen in the IRAC colors 
in addition to PAH line emission. 

4.2. 3.6 micron radial profile 

The optical morphology of NGC 6240 is highly irreg- 
ular. This is the product of two effects: uneven dust 
obscuration and actual morphological irregularity. Our 
8 /im observations confirm the clumpy nature of the dust 
(see Figure [3]) . With a reduction of the extinction by a 
factor of >10 over the optical, IRAC 3.6 and 4.5 /xm data 
is able to penetrate the dust obscuration. To this end, 
we use the elliptical isophotes discussed in § 14.11 to cre- 
ate a radial profile for NGC 6240 in 3.6 /xm to determine 
the true stellar distribution. We first fit an r^^^ law, 
shown in Figure [51 to determine how well this remnant 
approximates an elliptical galaxy . A similar calculation 
was done by iDovon et al.l (119941 ) in the K-band. Our 
results and theirs are very similar, but their K-band res- 
olution allows them to extend the profile inside of 1 kpc. 
IRAC's resolution does not allow this but its sensitivity 
allows us to extend the p rofile to 25 kpc, over twice as 
far as iDovon et al.l ()1994l ). Outside of 5 kpc, the fit ap- 
pears to agree well with a r^/'* law. The profile shows 
some "excess light," a deviation above an r^^^ law, in 
the central «1.5 kpc of NGC 6240. I n a s ample of 51 
merging galaxies, iRothberg fc JosephI (|2004l ) found that 
16 showed similar indications of excess light. They pro- 
pose three possible explanations. First, that this could be 
from a population of stars formed in a starburst triggered 
by dissipative collapse during the merge r as is shown in 
the models of iMihos fc Hernguis^ (|1994D . Secondly, that 
this light is due to a central AGN. Finally, this light could 
be due to the nuclear struct ure of the progenitor gal axies 
as shown by the models of iFulton fc Barnea (|2001h . In 
the case of NGC 6240, we know that two central AGN 
exist, making an AGN at l east part of t he explanation. 
By comparing our profile to lDovon et al.l (1994) it is clear 
that our profile starts outside the peak of either nuclei, 
so we cannot attribute this light to either nuclei. It is 
still unclear whether this light is from stars that were 
originally in one of the galaxies or part of a young star- 
burst populati o n. How ever using numerical simulations, 
iHopkins et al.l (|2008bl ). show that if extra light compo- 
nents are robustly separated out from the light profile 
of the galaxy, merger remnants ubiquitously contain a 
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Fig. 8. — The 3.6 /xm surface brightness profile derived from the 
elliptical isophotes plotted against r^'^. This fit gives an efi^ective 
radius of ~ 8 kpc. The dashed line is the effect of the PRF. 

component created primarily by the merger-driven star- 
burst, that is reflected in this extra light. Given this 
result, it is most likely that the "extra- light" seen in 
NGC 6240 is primarily from a starburst component re- 
sulting from the merger. This extra light component 
has some dependence on the progenitor galaxies initial 
gas fraction (though, counter-intuitively, the extra light 
decreases with increasing progenitor gas mass because 
at higher initial gas masses more gas is consumed in a 
starburst at first passage), bu t the relationship has high 
scatter for individual objects (iHopkins et al.ll2008bl ). In 
samples examined by IHopkins et al.l ( 2008bf ) this compo- 
nent's mass spans 3-20%. 

It seems surprising that an active merger is well fit 
by an r^'^ law, especially when its western half appears 
so disky. To check for signatures of a remnant disk, we 
averaged two radial cuts separated by one IRAC pixel 
along the major axis of the remnant. The locations of 
the cuts and the resulting profile are shown in Figure l9l 
On the northern side of the galaxy (crosses in Figure IHl), 
the profile appears exponential beyond 10 kpc. We fit 
this portion of the profil e with the edge-on gala x y sur - 
face brightness model of Ivan der Kruit fc Searlg ( 19811 ) 
(solid line in Figure[9]) which matches the data very well. 
The southern side is clearly not simply exponential. We 
attempt to fit the south profile beyond 5 kpc (the dashed 
line in Figure IH]). This fit only agrees well for the portion 
of the profile from 5 to 15 kpc. This is not surprising 
given the truncated and asymmetric morphology of the 
southern side of the remnant. This, as all other aspects 
of this remnant, reflect a galaxy in transition, in the pro- 
cess of relaxing while some portions of the remnant disk 
remain. 

5. CONCLUSIONS 

NGC 6240 is an object experiencing a rarely observed 
stage of evolution, showcasing complex combinations of 
processes. We add Spitzer-IRAC imaging and photome- 
try to the active study of this object. We show that stars 
and dust have very different distributions in the remnant, 
both regionally and radially from the central AGN. We 
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Fig. 9. — Left: The two major axis cuts averaged to make the profile in Figure [9] overlaid on the 3.6 ^m image with contours. Right: The 
surface brightness profile derived from the major axis cuts in the left panel. The red crosses are the northern side of the galaxy and the 
blue circles are th e southern side of the galaxy. The outer parts look exponential and are fit with a theoretical surface brightness profile 
llvan der Kruit fc^ Searle 1981). These fits are the solid red line on the north side and the dashed blue line on the south side. These indicate 
the remains of a disk in this remnant. The black dotted line is the 3.6 /xm PRF. 



examine how outflows in the X-ray, Ha and CO correlate 
with 8 /im PAH emission. We show that though, when 
radially averaged, the stellar profile follows an r^/'* law 
surprisingly well, regions of the galaxy still show disk- 
like stellar profiles. All this supports the basic picture 
of a galaxy experiencing the effects of a major merger, 
which treats stars and gas/dust very differently, and can 
create a gas poor spheroid from a gas rich spiral galax- 
ies. Though we can say little about the progenitors of 
this system, the sensitivity of IRAC data to low surface 
brightness emission adds detailed knowledge of the dis- 
tribution of stars and gas in this remnant. This provides 
constraints for large samples of hydrodynamic simula- 
tions coupled with radiative transfer which, if they are 
able to match the characteristics of NGC 6240, may be 
able to determine its history and future. 
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